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A Basis  for  Homogeneous  Polynomial  Solutions 
to  Homogeneous  Constant  Coefficient  PDE’s: 
An  Algorithmic  Approach  through  Apolarity 


Michel  Pocchiola  and  Gert  Vegter 


Abstract.  Some  recent  methods  of  Computer  Aided  Geometric  Design 
are  related  to  the  apolar  bilinear  form,  an  inner  product  on  the  space  of 
homogeneous  multivariate  polynomials  of  a fixed  degree,  already  known  in 
19th  century  invariant  theory.  A generalized  version  of  this  inner  product 
was  introduced  in  [8]  to  derive  in  a straightforward  way  some  of  the  recent 
results  in  CAGD.  Here  we  extend  this  work  by  applying  it  to  compute 
solution  spaces  of  homogeneous  constant  coefficient  PDE’s. 


§1.  The  Homogeneous  Apolar  Bilinear  Form 
1.1.  Review 

In  [8]  we  introduced  a generalization  of  the  apolar  bilinear  form  defined  on 
the  space  of  homogeneous  polynomials  (of  a certain  degree,  and  with  a fixed 
number  of  variables).  This  bilinear  form,  used  extensively  in  the  symbolic 
method  of  the  classical  theory  of  invariants,  has  been  revitalized  by  Rota  and 
his  co-workers,  cf  [2]  and  [4].  In  CAGD,  a similar  binary  form  on  the  space  of 
univariate  polynomials  of  a fixed  degree  has  been  studied  by  Goldman  [3].  It 
is  related  to  the  blossoming  approach  introduced  by  Ramshaw  [7]. 

In  this  section  we  review  some  of  the  properties  of  the  apolar  bilinear 
form.  Then  we  extend  [8]  by  studying  constant  coefficient  partial  differential 
equations  of  the  form  p{d)f  = 0,  where  p is  a fixed  multivariate  homoge- 
neous polynomial.  In  particular,  we  derive  an  algorithm  computing  a basis 
for  solution  spaces  consisting  of  homogeneous  polynomials  of  a fixed  degree. 
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1.2  Vector  spaces  of  forms 

Let  ei,. ..  ,es  be  the  standard  basis  vectors  on  ]RS,  and  let  x = (xr, . . . ,xs ) 
be  the  standard  coordinates  on  1RS.  The  standard  inner  product  on  1RS  is 
denoted  by  (•,  •),  i.e.,  (u,  v)  = ulVl  +■■■  + usvs,  for  u,  v £ IRA 

A central  object  in  this  paper  is  the  space  of  real  homogeneous  poly- 
nomials of  degree  n on  Rs,  denoted  by  7fn(Rs).  A polynomial  in  7fn(lRs) 
is  the  sum  of  monomials  of  the  form  cax“'  where  cQ  € 1R  and  a = 

(ai, . . . , as)  £ is  a multi-index  of  weight  |a|  = «!  + •••  + as.  For  con- 
venience the  monomial  sc"1  ■ • • j;“’  is  denoted  by  x° . Linear  homogeneous 
polynomials  on  IRS  are  of  the  form  f(x)  = {u,x),  for  some  u £ IRA  We  denote 
/ by  («,•)• 

For  multi-indices  a = (07, ...  ,as)  and  (3  = ((3U...  ,/?s)  in  ZZ|0  we  define 
a | /?  iff  ai  < ^ for  i = 1,  • • • , s.  The  relation  | is  a partial  order  on 
Note  that  a \ f3  iff  there  is  a A € 2Z>0  such  that  /?  = q + A.  We  shall  write 
a-\-p  if  it  is  not  the  case  that  a | (3.  A monomial  order  is  a linear  ordering 
<mon  on  2Z|0  such  that  (i)  if  a<mon/3  and  7 £ 2Z>0,  then  a + 7<mon/3  + 7, 
and  (ii)  <mon  is  a well-ordering  on  K>0,  i.e.,  every  non-empty  subset  of  2Z>0 
has  a smallest  element  with  respect  to  <m0n-  We  use  the  notation  a<mo^/3 
in  case  a<monP  and  a ^ fj.  Furthermore,  we  use  the  property  that  a<mon6 
whenever  a \ /3.  Well  known  examples  are  the  graded  (reverse)  lexicographic 
orders,  defined  by  a<monp  if  |a|  < \0\,  or  |q|  = \f3\  and  in  a - (3  the  left- 
most (right-most)  non-zero  entry  is  negative  (positive).  Monomial  orders  play 
a paramount  role  in  algorithms  for  multivariate  polynomials,  especially  with 
regard  to  termination  conditions;  See  e.g.  [1] . 

The  set  of  multi-indices  in  ZZ>0  of  weight  n,  denoted  by  r5i„,  is  a finite 
set  with  #rs,n  = (n+*  J)  elements.  For  a £ rSin  the  factorial  function  is 
defined  by  a!  = 07!  • ■ • as!,  and  the  multinomial  coefficient  (")  is  defined  by 


With  a polynomial  f(x)  = J2aers  m c<* xa > we  associate  the  homogeneous 
differential  operator  f(d ) = Eaers  m ca  0° , where  da  = d*'  Here 

d = (di, . . . ,ds),  with  di  = d/dxi . The  directional  derivative  Du  : ?fn(]Rs) 
^n-i(IRs)  with  respect  to  u £ 1RS  is  the  differential  operator  (u,d),  i.e., 
Du  — ui  d\  + ■ ■ ■ + us  ds.  Note  that  di  = (e^  d)  = De< . Considering  e*  as  a 
multi-index  of  weight  one,  we  also  have  di  = dei . 

1.3.  Apolar  pairing 

This  subsection  is  concerned  with  a straightforward  generalization  of  the 
rather  well-known  apolar  inner  product  [/,  5]  = f{d)g , defined  on  the  space 
of  homogeneous  polynomials  TinflR’).  The  main  result  concerns  a character- 
ization of  this  inner  product  in  terms  of  three  simple  properties  that  will  be 
the  basis  for  the  construction  of  special  bases  of  7fn(lRs)  in  later  sections. 
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Definition  1.  For  fixed  integers  m and  n,  with  0 < m < n,  the  apolar  pairing 
is  the  map 

[ • , • : Wm(ffts)  x Wn(Ks)  - H„-m(lRs), 

associating  to  the  homogeneous  polynomials  f € 7fm(Rs)  and  g £ 7in(Rs) 
the  homogeneous  polynomial  [f , g]m,n  of  degree  n — m,  defined  by 

(n  - m)!  ... 

[/ ! s]m,n  — f(d)g- 

Note  that  we  have  in  fact  a family  of  pairings,  one  for  each  pair  of  integers 
m and  n with  0 < m < n.  In  this  paper,  the  term  pairing  refers  to  the  whole 
family  of  bilinear  maps.  From  now  on  we  shall  drop  the  subscripts  m and  n, 
since  they  are  implicitly  known  as  the  degree  of  the  first  and  second  argument 
of  the  pairing  operator. 

Theorem  2.  The  apolar  pairing  is  the  unique  bilinear  pairing  with  the  fol- 
lowing properties: 

1)  (Apolar  pairing  with  constants).  For  f £ 7f„(lRs): 

[1,/]  = /, 

where  1 € 7fo(Rs)  is  the  constant  homogeneous  polynomial  of  degree  0. 

2)  (Apolar  pairing  with  linear  forms).  For  f £ 7f„(Its)  and  u € Hs: 

[(u,.),f]  = ±Duf. 

3)  (Transposition  of  a homogeneous  factor).  For  f\  € ?fmi(Rs),  fi  £ 
7fra2(lRs),  and  g £ Ttn{Rs),  with  mi  + m2  < n: 

[/1/2 ! 5]  = [/1 ) [h  > 9]  ]• 


It  is  obvious  that  apolar  pairing  is  a bilinear  operator,  satisfying  these 
properties.  For  the  proof  of  uniqueness,  we  refer  to  [8].  Identifying  the  space 
of  zero  degree  polynomials  with  1R,  we  see  that  for  n = m,  apolar  pairing 
corresponds  to  a real  bilinear  form  on  the  space  of  homogeneous  polynomials 
of  degree  m.  The  next  result  states  that  this  bilinear  form  is  even  an  inner 
product.  Again,  the  proof  is  contained  in  [8]. 

Proposition  3.  The  apolar  bilinear  form  [ • , ■ ] : 7dm(lR's)  x 7fm(]R5)  — * 1R  is 
an  inner  product  on  the  space  of  homogeneous  polynomials  of  degree  m. 
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1.3.  Dual  bases 

First  we  recall  the  definition  of  a dual  basis  pair  with  respect  to  the  apolar 
inner  product  [■,  •]  on  ?tm(Rs). 

Definition  4.  The  dual  basis  of  a basis  {fa  \ a £ r9m}  of  7im(Rs)  is  a 
collection  {ga  \ a £ rsm}  of  polynomials  in  ?fm(] Rs)  such  that,  for  a,/3  £ 

r 

1 s3m* 

[ fa  i 90  ] = 

It  is  easy  to  prove  the  standard  fact  from  linear  algebra  that  a dual  basis 
is  indeed  a basis.  Given  a dual  basis  pair,  a polynomial  / £ 7im( Rs)  can  be 
expressed  with  respect  to  either  basis  in  terms  of  coefficients  depending  on 
the  other  one: 


/=  [ 9a  i / }fa  ~ [/«’/]&*•  (x) 

oers,m 

Example.  (Dual  of  homogeneous  Bernstein-Bezier  basis).  Let  {x1,...,®8}  be 
a basis  of  Rs,  and  denote  by  «i(x),  ...,u3(x)  the  coordinates  of  any  x £ Rs 
in  this  basis.  The  polynomials 


where  a £ rs,„,  form  the  homogeneous  Bernstein-Bezier  basis  of  7fn(Rs) 
with  respect  to  the  basis  {x1,...,®5}  of  Rs.  Its  dual  basis  consists  of  the 
polynomials 

b(v)  = (xl>y)01  •••(*’,  j i)0’, 

i.e.,  [Ba  , Ip]  = ba,p.  For  the  proof,  see  [8]. 

§2.  Solving  Homogeneous  Constant  Coefficient  PDE’s 

We  now  show  how  dual  bases  can  be  used  for  the  efficient  computation  of 
a basis  for  the  solution  space  of  a homogeneous  partial  differential  equation 
with  constant  coefficients,  i.e.,  the  space 

{/  £ Hn( Rs)  | p(d)f  = 0}.  (2) 

Here  p £ ?fm(Rs)  is  a polynomial,  p / 0,  that  will  be  fixed  throughout 
the  paper.  Furthermore,  m and  n denote  fixed  integers  such  that  0 < m < 
n.  Our  approach  is  both  an  alternative  and  an  algorithmic  counterpart  of 
Pedersen’s  work  [5,6].  These  papers  deal  with  algebraic  properties  of  the  space 
of  solutions.  We  continue  Pedersen’s  work  by  presenting  an  optimal  algorithm 
for  the  computation  of  a basis  for  the  solution  space.  Our  techniques  are  new, 
since  they  are  based  on  properties  of  dual  bases,  together  with  some  recursive 
properties  of  the  apolar  bilinear  form  introduced  in  [8]. 


329 


Apolar  Bilinear  Form  and  PDE’s 


2.1.  Characterizing  a basis  for  the  space  of  solutions 

From  now  on  we  consider  a family  of  functions  {/<*  | a G Z>0}  such  that 
(i)  fa'fg—  fa+g i and  (ii)  for  all  n > 0,  the  set  {fa  | a G FSi„}  is  a basis 
of  7fn(lRs).  The  dual  basis  of  the  latter  set  is  denoted  by  {ga  \ a G rs>n}. 
An  example  of  such  a pair  of  bases  is  formed  by  the  Bernstein-Bezier  basis, 
together  with  the  lineal  polynomials  introduced  in  the  example  at  the  end  of 
the  preceding  section. 

Lemma  5.  For  a, (3  G 7Zl,0  with  |o|  < \j3\: 


[fa  i 90 


90— a i 

o, 


if  a | (3, 
if  a-j-/3. 


Proof:  Let  m = |a|  and  n = |/?|,  and  let  / — [fa,gp]  € 7f„_m(lRs). 
Consider  the  apolar  inner  product  [/7 , /],  for  7 6 rs>n_m.  Since  fa  ■ f1  = 
fa+y,  transposition  of  a factor  (See  Theorem  2,  part  3)  /„  yields:  [/7 , /]  = 
[fa+y , 9p } = 6a+y,0-  First  consider  the  case  a-\-(3.  Then  and 

hence  [/7  , / ] = 0,  for  all  7 G rSi„_m.  Since  the  apolar  pairing  is  an  inner 
product  and  {/7  | 7 G FSi„_m}  is  a basis  of  it  follows  that  / = 0 

in  this  case.  If  a \ f3  the  previous  derivation  shows  that  [/7 , /]  = 67 ,/?-<*,  so 
identity  (1)  implies  f = E7era  [f,fy]9y  = 9p-a-  □ 

In  the  following,  our  fixed  polynomial  p in  (2)  is  given  in  the  form 


P=  Ca  £*>  where  C<*  = \V,  9a)- 

aer3,m 


The  following  result  characterizing  the  kernel  of  a polynomial  differential  op- 
erator is  the  key  ingredient  for  the  algorithm  developed  in  the  next  section. 

With  p we  associate  the  linear  map  Dp  : 7f„(]Rs)  — > TLn-m(Mf)  defined 
by  Dp(f)  = [p,  /],  and  the  map  Tp  : ?fn_m(] Rs)  — > 7fn(Rs)  is  multiplication 
by  p,  i.e.,  Tp(f)  = p ■ f . Given  an  integer  k,  and  a subspace  U C 7-fjt (3Fbs), 
we  denote  by  t/x  the  orthogonal  complement  of  U with  respect  to  the  apolar 
inner  product  [• , •]  on  Hk(Rs)- 

Proposition  6. 

1)  KerDp  = ( ImTP)x . 

2)  The  map  Dp  is  onto. 

Proof:  Theorem  2,  part  3,  implies  that  Tp  and  Dp  are  adjoint  operators, 
i.e.,  [Tp(f) , g]  = [f , Dp(g)],  for  / € Hn_m(lRs)  and  g g 7f„(lRs).  The  first 
claim  follows  from  this  identity.  Now  since  Tp  is  injective,  the  result  of  the 
first  part  implies  that  dimKerDp  — dimTtn(TRa)  — dim7fn_m(lRs).  Therefore, 
dimlmDp  = dim ?tn_m(IRs),  and  hence  Dp  is  onto.  □ 

As  a special  case,  consider  the  polynomial  p = fao  for  some  ao  € rs?m. 
According  to  Lemma  5,  Ker Dp  contains  gp  whenever  (3  G rSi„  such  that 
OiQ-\-[3.  Since  oq  I (3  iff  (3  is  of  the  form  /3  = oq  + A for  some  A g rs>n_m,  it 


330 


M.  Pocchiola  and  G.  Vegter 


follows  that  #{/?  G r3,n  | ao-\-(3}  — #rs,n  - #r s,n-m  = dimKerDp.  The  last 
equality  follows  from  Proposition  6,  part  1.  Therefore,  a basis  for  the  solution 
space  Ker Dp  is  the  collection  {gp  \ f3  € rs  n and  op  -(-/?}•  The  following  result 
generalizes  this  special  case. 

Theorem  7.  (Basis  for  solution  space  of  PDE).  Let  <mon  be  a monomial  order 
on  E>0  and  a0  € rs,m  be  defined  by  a0  = min<mon{a  G TSi7n  | [p,  ga]  ^ 0}. 
Furthermore,  for  any  A € TS]n_m,  let  p\  G Wn-m{  Rs)  be  the  polynomial 
defined  by  px  = [p,  gao+x}-  Then 

1)  The  set  Vn-m  = {p\  \ X G rSin_m}  is  a basis  of  7fn_m(lRs). 

2)  Let  Qn-m  = {<7a  1 A G rs,„_m}  be  the  dual  of  the  basis  Pn-m  of 

Rs),  i.e.,  [pa  , qp  ] = 6Xp.  A basis  for  the  solution  space  KerDp  = 
{/  G 7in{\ R5)  | p(d)f  = 0}  is  the  set 

{§p  I P € Ts,„  with  q0+/3},  (3) 

where  cjp  G Hn{  Rs)  is  defined  by 

90  = 9p~  X]  [P  ’ 9a  , ff/3  ] ,9a0+A- 

^€r3>n  — m 

Remark.  The  first  claim  of  Theorem  7 is  not  necessarily  true  for  other  choices 
of  a o-  Consider  e.g.  the  polynomial  p(x)  = 2xf  + 2aqx2  + x2,  and  let  ao  — 
(1, 1)  G T2,2-  Here  we  take  the  monomial  basis  for  the  space  of  polynomials  of 
degree  n on  R2,  i.e.,  we  take  fp{x)  = , for  x G R2  and  |^|  = /3i  + /?2  = n. 

The  dual  basis  consists  of  the  functions  gp,  where  gp{x)  — Q)xfl.  For  A G T2,2 

we  have  gao+x(x)  = (a^A)xj1+1X22+1.  Take  q(x)  = 2x2  — 2xjX2  + x2,  then 
p(x)  ■ q(x ) = 4xi  + X2,  and  hence,  for  all  A G 1^2,2: 

[Pa  , q]  = [[p,  Pao+A  ] , q]  = [p-q,  Poo+a]  = 0, 

yet  q / 0.  Hence  the  functions  px,  where  A ranges  over  T 2,21  do  not  constitute 
a basis  for 

To  prove  Theorem  7 we  need  the  following  two  lemmas. 

Lemma  8.  For  A,  p € rSi„_m  we  have 

r»,  / 1 _ / cao+A-/ii  JTp  | ao  + A, 

1 /r*  I |0,  otherwise. 

Proof:  The  proof  consists  of  a straightforward  calculation: 

[PA,  U]  = [fti,  [p,  Pa0+A  ] ] = [p,  [U,ga  o+A  ] ] 

_ f[P,  Pao+A-^],  if  P | QO  + A, 

1 0,  otherwise. 


The  last  identity  is  justified  by  applying  Lemma  5.  □ 
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Lemma  9.  Let  Ao  € rs,n-m.  For  f € ?fn_m(Rs),  the  following  statements 
are  equivalent: 

1)  [pa  , /]  = 0,  for  all  A<  mon  Ao, 

2)  [Pa  1 /]  = o,  for  all  X<monX0. 

Proof:  It  follows  from  (1)  that 

[pa./]=  ^2  [PA.  U][9n, /]•  (4) 

ner, 

,n  — m 

Consider  p £ rS)„_m  such  that  A<mon/i,  then  ao  + A — /i<monao-  Therefore, 
the  definition  of  ao  implies  cao+A-/i  = 0.  In  view  of  Lemma  8,  we  know  that 
[px  , fji  ] is  equal  either  to  cao+\-^  or  t°  0»  so  in  any  case  we  have  [pa  , } =0. 

This  observation  allows  us  to  write  (4)  as 

[pa  , /]  = cao[gx  , f]  + [pa>  U][9ii,  /]•  (5) 

This  identity  shows  that  the  first  statement  implies  the  second  one.  So  assume 
that  statement  2 holds.  We  may  assume  that  / / 0,  otherwise  there  is  nothing 
to  prove.  Let  Ai  be  the  least  multi-index  with  respect  to  the  monomial  order 
<mon  such  that  [gXl  , f]  + 0.  Then  (5)  implies  [pAl  , /]  = cao[gXl  , f } is 
nonzero.  Hence  Ao<m0nAi.  Consequently  [g\  , f]  — 0 for  all  A<monAo,  which 
is  statement  1.  □ 

Proof  of  Theorem  7:  Let  U C Wn_m(]Rs)  be  the  space  spanned  by  the 
Pa,  A € Ts,n_m.  Since  = dimH„-m(l Rs),  it  is  sufficient  to  prove 

that  U = 7fn_m(R3),  or,  equivalently,  that  U±  = {0}.  Thus,  if  / € UL , 
then  [/,  pa]  = 0,  for  all  A € TSi„_m.  According  to  Lemma  9 this  implies 
[/,  pa]  = 0,  for  all  A € TSin_m,  so  / = 0.  This  proves  1).  Now  in  view 
of  Proposition  6,  the  space  Ker Dp  is  of  dimension  #TSi„  — #TSi„_to,  i.e.,  of 
dimension  #{/3  € r5]„  | ao -]-/?}.  On  the  other  hand,  it  is  straightforward  to 
see  that  the  polynomials  p^,  (3  € TS:„  with  ao-f/3,  are  linearly  independent. 
Therefore,  in  order  to  prove  that  they  form  a basis  of  KerDp  we  just  have  to 
prove  that  they  belong  to  KerDp.  Taking  Proposition  6,  part  1,  into  account, 
we  actually  have  to  check  that  [p  ■ <7P  , Pp  ] = 0,  for  all  /3  £ rs,n  with  ao-f/3. 
Since 


7ers,„ 

= Y2  [[p-  5«0+a],  PM]/ao+A  + J2  [p-Pm>  9~,]fy 

ao+7 

= fct0+n  + yi  [P  ' i ] I'll 
7er,,„ 

ao+7 

it  follows  that  [p  ■ pM  , pp ] = 0,  for  all  (3  € TS:n  with  ao-f/3.  □ 
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2.2.  Computing  a basis  for  the  space  of  solutions 

We  now  present  a simple,  efficient  algorithm  for  computing  the  dual  basis 
Qn-m,  as  well  as  an  example  showing  how  the  algorithm  works.  Recall  that 
for  a = (ai, . . . ,os)  G r«,mi  the  number  ca  is  equal  to  [p,  ga  ].  We  extend 
this  definition  to  a G 7LS  by  putting  ca  = 0 in  case  at  least  one  of  the  entries 
ai,...,as  is  negative. 


Corollary  10. 

1)  The  dual  basis  Qn_m  = {qM  j p G rS]„_m}  ofVn—m  is  defined  recursively 

by 


Qp  — (//i  'y  Cao+v-iiQv)  ■ 

Cao  ver„„_m 


2)  For  fi  G rs,n,  with  oq  -f/3,  the  basis  function  G 7fn(Rs),  is  of  the  form 


90  — 90  y \ aft0  9a0+ti‘ 

where  the  coefficients  a^p  are  defined  recursively,  for  p G by 


cao+v~tiav0)  ■ 

a,n  — m 
M^mon  ^ 


Proof:  Recall  that  we  are  looking  for  a set  of  functions  Qn-m  = { q ti  | p € 
rs,n-m}i  such  that  [p„  , <jp]  = <5„M.  In  particular,  according  to  Lemma  9 the 
functions  satisfy  [gv  , q^}  = 0,  for  v<mOn0-  Therefore,  G Span{/„  | v 6 
r s,n-m  and  p<m0n"},  or,  equivalently: 


q n € Span({<p,  | v G Ts,n_m  and  \i<monv}  U {/,,}).  (6) 


Assume  we  have  determined  qv  for  v G rs,n_m  with  p<m0 n^-  To  compute 
q M satisfying  (6),  we  have  to  determine  constants  d^v,  for  p,^  G TSi„_m  with 
M<mon^!  such  that 

q^  — dppfp  + y ^ dpxqx- 

A£T 5 ,n  — m 
M^mon  ^ 

Since  Vn-m  and  Qn-m  are  dual  bases,  the  constants  dliu  are  uniquely  deter- 
mined by  the  condition  [p„  , Combining  the  last  two  identities  we 

see  that 

[Pi/ , ] = d^[p„  , ] + y d^x6x,v 

3,n-m 
^^mon  A 

From  this  identity,  which  holds  for  all  v G rSin_m  with  p<mo nv,  we  derive 

*W  = - , 

ca  o 

J CotQ  — ll  + l/  c , 

d^t/  — , IOr 

Ca0 
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which  proves  the  first  part.  Now  put  a^p  = [p-gM  . gp]  in  (3).  Then,  according 
to  part  1, 


av(3  ~ (I?'  fn>9p]  ^ 1 ca0+v-l^[p  ' Qv  > 9p  ])  • 

*er3,„_m 

fJ-  ^ mon  V 

Since  [p  ■ /M , gp  ] = \p , [/p  , gp  ] ] = [p , gp-^  ],  the  proof  is  complete.  □ 

The  algorithm  for  computing  a basis  for  the  solution  space  of  the  partial 
differential  equation  p(d)f  = 0 is  now  simple: 


Algorithm  (for  computing  a basis  for  Ker£>p). 


forall  p € r\n_m  in  decreasing  <mon-order  do 
forall  0 eTsn  with  ctQ-\-0  do 


^ ' cao +u—y.avp\ 


Example.  Consider  on  1R3  the  homogeneous  constant  coefficient  PDE 

d2f  d2f  82f 

8x38x2  8x\  6x28x3  ’ 

corresponding  to  the  homogeneous  polynomial  p(x)  = X1X2  — x%  + X2X3.  In 
particular,  the  setting  of  this  example  corresponds  to  s = 3 and  m = 2. 
We  determine  a basis  for  the  solution  space  in  /^(IR.3),  i.e.,  we  take  n = 3. 
To  this  end,  consider  the  graded  reverse  lexicographical  order  on  Z>0.  Let 
fp{ x)  = x13,  and  let  gp(x)  = (fyx13,  where  n = \0\.  In  this  example  we  denote 
functions  indexed  by  a = (i,j,  k)  € 0,  like  fa,  by  /y*,. 

The  sets  { fp  \ 0 € TSi„}  and  {gp  | 0 € rs>n}  are  dual  bases,  and  moreover 
fa  • fp  = fa+p,  so  the  conditions  for  applying  Theorem  7 and  the  algorithm 
from  this  section  are  satisfied.  In  the  notation  of  Theorem  7,  we  have  ao  — 
(1,1,0).  Note  that  p = fno  - /020  + /oil,  so  cn0  = 1,  c02o  = -1  and 
coil  = 1,  whereas  all  other  coefficients  Cy*,,  with  i +j  + k = 2,  are  zero.  Now 
the  coefficients  a^p  are  computed  according  to  the  algorithm  above,  in  other 
words  we  successively  determine  the  rows  in  the  following  table  (computing 
for  each  row  the  entries  in  arbitrary  order): 


anP 

0 = 300 

201 

102 

030 

021 

012 

003 

p — 001 

0 

0 

0 

0 

-1 

1 

0 

010 

0 

0 

0 

-1 

1 

0 

0 

100 

0 

0 

0 

-1 

2 

-1 

0 

This  table  corresponds  to  the  following  seven  basis  functions  of  Ker Dp:  <7300, 
3201,  3102,  3030 +3120 + 3210,  302i+3in  -3120  — 23210,  3012  - 3in +3210,  and 
3003- 
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These  functions  can  be  turned  into  monomial  form  by  straightforward  substi- 
tion,  yielding  the  following  explicit  basis  for  the  solution  subspace  of  7/3  (R3): 
x3,  3xj%3,  3x10:3,  3xf*2  + 3xi*2  + X3,  —6XjX2  — 3X1*2  + 6x1X20^3  + 3x3X3, 
3xfx2  — 6X1X2X3  + 3*2X3,  x3. 


Acknowledgments.  We  are  greatful  to  the  anonymous  referee  for  corrections 
and  suggestions  for  improvement. 
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